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Abstract 

A set of one-loop vertex and box tensor-integrals with massless internal particles has 
been obtained directly without any reduction method to scalar-integrals. The results with 
one or two massive external lines for the vertex integral and zero or one massive external 
lines for the box integral are shown in this report. Dimensional regularization is employed 
to treat any soft and collinear (IR) divergence. A series expansion of tensor-integrals with 
respect to an extra space-time dimension for the dimensional regularization is also given. 
The results are expressed by very short formulas in a suitable manner for a numerical 
calculation. 
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1 Introduction 



The LHC (Large Hadron Collider) project 1 at CERN is planned as the next-generation energy- 
frontier experiment. One of the physics motivations of LHC experiments is to discover the Higgs 
particle, which is the only one missing ingredient in the standard model. In order to establish 
the model completely, it is essential to find it and to investigate in its details. It should also be 
important to search for new phenomena beyond the standard model through any tiny deviation 
in experimental observations from theoretical predictions. Further, not only searching for new 
phenomena, but also performing precise measurements of parameters included in the standard 
model is another important issue of LHC. 

LHC has employed colliding proton-proton beams in order to achieve beam energies as high 
as possible, which should enhance the possibility to find new particles/phenomena. However, 
a proton machine must have a large QCD background, since the proton is a composite particle 
constructed by strongly interacting particles, such as quarks and gluons. This veils signals 
with large backgrounds. In order to extract as much physics information as possible from 
experimental data contaminated by huge backgrounds, the behavior of the background should 
be understood in detail. This implies that one should precisely understand QCD, because it 
entirely governs the background. 

The background of any proton-proton colliding experiment must be completely predicted 
by QCD, in principle. However, in fact, it is very a difficult task to make precise predictions 
because of the large coupling constant of QCD. Moreover the lowest level (tree-level) calculation 
does not have any predictive power for the event rate, because there is no good renormalization 
point well-defined experimentally. We need higher order perturbation calculations for precise 
predictions of the background behavior. 

Loop integration is one of the critical issues for computing these higher order corrections. 
In general, N-point tensor- and scalar-integrals with massless internal lines including infrared 
(IR) divergence must be calculated in QCD. Since an arbitrary number of dimensions must 
be used in QCD to regulate any IR divergence, the usual method for the standard model 2 
cannot be used directly. Various methods to reduce (N > 5)-point integrals into (N — l)-point 
integrals .4_ with a dimensionally regulated scheme are proposed. Then, all of the necessary 
(N > 5)-point integrals can be expressed by a linear combination of box (four-point) tensor- 
integrals. Usually, box tensor-integrals are further reduced to four or less point scalar-integrals, 
and then numerically evaluated to obtain higher order corrections. The IR finite box integrals 
are obtained in [JJj; for the IR divergent case, box integrals with zero or one external mass 
are given in 0|. All IR divergent box integrals are treated in using the partial differential 
equation method. Another approach to all possible box scalar-integrals with massless internal 
lines with zero- to four-external massive lines is proposed for the IR divergent casejH] and the 
IR finite case|0.. Recently, two independent formalisms were proposed for calculating one-loop 
virtual corrections with an arbitrary number of external legs jlUl ITT] . 

We propose a new method to calculate tensor-integrals directly, and not to use a reduction 
method from tensor-integrals to scalar ones in this report. A set of one-loop vertex and box 
tensor-integrals with massless internal particles is given in terms of hypergeometric functions. 
Dimensional regularization is employed to treat any IR-divergence. A series expansion of tensor- 
integrals with respect to an extra space-time dimension in the dimensional regularization is also 
given in this report. All results are expressed by very short formulas with a suitable manner 
for numerical calculations. 

A general form of three-point tensor-integrals is given in section 2, and that of four-point 
ones in section 3. The obtained results are numerically checked in section 4. A numerical 
calculation method for the hypergeometric function and their series expansion with respect to 
an extra space-time dimension is given in Appendix A. The series expansion of the general 
form of tensor-integrals can be represented in terms given in Appendix A. Those results are 
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given in Appendix B. 



2 Vertex integral 

2.1 Massless one-loop vertex integral in n-dimension 

The tensor-integral of a massless one-loop vertex with rank M < 3 in a space-time dimension 
of n = 4 — 2euv can be written as 

\ £ uv f d n k k^ • • • k v 




(/4) £ 



(27r) n D1D2D3' 



where 



D x = A; 2 + iO, 

£>2 = (A: + p 2 ) 2 + i0, 

#3 = (fc+P2+P3) 2 +^0, 

and pi is the four momentum of an z'th external particle (incoming), k^ the loop momentum, 
and \xr the renormalization energy scale. An infinitesimal imaginary part (iO) is included 

(3) 

to obtain analyticity of the integral T^... v . Momentum integration can be performed using 
Feynman's parameterization, which combines propagators. After momentum integration, an 
ultra-violet pole is subtracted under some renormalization scheme. Then space-time dimension 
is replaced as ejjv ~~ * ~~ £ IR to regulate an infrared pole. Finally, the tensor-integral is expressed 
in the following form|12j: 

T (3) r r i jifjl 2 2 (t) (i)s 



4iplplph^,rS) = r^Tdy^y^ 



where 



(4vr) 2 (47r^) e « 7o Jo D 1 '^ ' 

D = fax - p 2 y) 2 - pxy - p\x - p\y - iO, 
P = Pl-{P1+P2?. (2.1) 

The masses of internal particles are assumed to be massless. The remaining task is to perform 
the parametric integration of J3. 

2.2 Two on-shell, one off-shell external legs 

For the case of two on-shell and one off-shell external particles, we set p\ = p\ = 0,p 2 7^ 
without any loss of generality. The integration can be done in a straightforward way: 

J 3 {0,0,pl;n x ,n y ) = 77^ ^77^ f/^ [ dx [ dy " 



(4tt) 2 (4vr^|) £ « J J (-p 2 xy - iO) 1 "^ 

e IR T(-e IR ) f-p3 2 Y IR 1 B(n x + e IR ,n y + sir) 



(2 2) 

(4tt) 2 \^P 2 r ) ~p\ n x + n y + 2e IR 

where P3 2 = p§ + iO, and B(-, •) is a Beta function. The infrared structure of the tensor-integral 
can be obtained by expanding Eq. (|2.2|) with respect to sir. The results of expansions under 
the MS scheme are shown in Appendix B. When both n x and n y are non-zero, there is no 
IR-divergence as 

9 n 1 (n r — — 1)! 

73(0,0,rf;„„„„) ' in j; n " y) , <««-<>). 
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2.3 One on-shell, two off-shell external legs 

In the case of one on-shell and two off-shell external particles, we set p\ = 0,p| 7^ 0,p| 7^ 0. 
The integral J3 becomes 

J 3 {Q,pl,pl;n x ,n y ) 

1 e IR r(-e m ) [ l dx f 1 ~ x dy _ * nx y ny 



(An) 2 (4vr^) £ « J J ({p 2 - p 2 2 )xy - p\ y(l - y) - iO) 1 -^ 

e IR T(-e IR ) ( -p 3 2 \ £lR 1 B(n x + e IR , n y + e IR ) 



(47r) 2 \4tth 2 r J -p\ n x + n y + 2e IR 

TP (\ 1 e 9 I - . jjzij ^ w *+ g ^ 
2^1 1, 1 - £/i?, 2 + n x ; — — — , 

V P3 / na; + l 

J 3 (0Aphn x ,n y ) g n (^^\ , (2.3) 



where 



0«O) = n + * IR 2 F 1 (l,l-s IR ,2 + n;z), (2.4) 
n + 1 

and 2-^1 (•> ■) ") is the hypergeometric function. A definition and some properties of the hyper- 
geometric function and its numerical evaluation can be found in Appendix A. When p\ — > 0, 
the hypergeometric function becomes 

n 4- 1 

2 Fi(l,l -sir, 2 + n x ;z) - — ^ (z -» 1). 

Then, the result Ea. ()2.3|) agrees with Ea. l|2.2jl when p\ — > 0. The infrared structure of the 
tensor-integral can be obtained by expanding Eq. (|2.3|) with respect to ej R . The results of 
expansions under the MS scheme are shown in Appendix B. 



3 Box integral 

3.1 Massless one-loop box integral in n-dimension 

Box integrations can be treated the same as in the vertex case. The tensor-integral of a massless 
one-loop vertex with rank M < 4 in space-time dimensions of n = 4 — 2ejjv can be written as 

„(4) _ / 2 \ £ uv f d n k fc^j ■ ■ k v 
L^Z ~ ™ J {2k)HD 1 D 2 D z D^ 

M 

where 

Di = k 2 + iO, 

D 2 = (k+ Pl ) 2 + iO, 

D 3 = (k + pi + p 2 f + iO, 

D 4 = (k + pi + p 2 + p 3 ) 2 + »0. 

After the same procedure as that used in vertex integration, we are left with following para- 
metric integration: 

r(2 — £l R ) f 1 r 1 ~ x rX—x—y x n X yn Vz n z 

Ji(s,t;plplplpl;n x ,n y ,n z ) = ^ J q dx ^ dy ^ 

(3.1) 
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where 



D = —sxz — t y(l — x — y — z) — p\ xy — p\ yz 

—p\ z(l — x — y — z) — p\ x(l — x — y — z) — iO, 
s = (pi+p 2 ) 2 , 
t = (pi+Pa) 2 - 

3.2 Four on-shell external legs 

When all external particles are on-shell (massless), p\ = p\ = p\ = p\ = 0, the integral of 
Ea.ipTT|l can be 

J4,(s,t;0,0,0,0;n x ,n y ,n z ) = 

!\ iyi , ^'•^ . (3.2) 



(4tt) 2 (4irfi R ) £lR J (-xzs - y(l - x - y - z)t - iO) 2 - £lR 

After applying the transformation 

x = r v, 

y = w (1 - r ), 

z = (1 — r) (1 — w), 

the integral becomes 

J 4 ( S ,t;0,0,0,0;n x ,n y ,n z ) = - \^ ~ £ f ] f dr r -l+«*+««(l - r )-i+n v +n z+em 

(4vr) 2 (4vr/i|) IR J 

dv I dw- 



2SIR 



o Jo (— s v{\ — w) — t (1 — v)w — iO) 

" W (4irfj, 2 R ) 6IR B(n x + e IR , n y + n z + e IR ) 



(4vr) 
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-l rl w n v(l-w) n *v n ° 
x I dv I dw- 



o J o 



(— s v(l — w) — t (1 — v)w — iO) 



The r-integral just gives the Beta function. Then, the u-integral can be done as 
I v = dv- 



eiR w sir-2 



t 2 (l + n x ) 



■2-eiR 



Fx 2 - e IR , 1 + n x , 2 + n x , -£ v 



o [(— s + sw + tw)v — tw — iO] 



where 



~ us 

Qw = J + ~ i 
t t W 

s = s + iO, 
i = t + iO, 

u = u + iO = (pi +P3) 2 + iO. 

Here, we use s + t + u = J2i Pi = 0- 

Further integration of the hype-geometric function is not straightforward. When a (or b) 
in the hypergeometric series 2-^1 (o> b, c; z) is a negative integer, the hyper geometric series is 
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truncated at some point, and becomes a polynomial. In oder to express our target integrand 
as a polynomial, a transformation formula, 



2.F1 (a, b, c; z) 



r(c)T(c-a-b) __ c 



T(c-a)T(c-b)' 



z a 2-F1 [a, a - c + 1, a + b - c + 1; 1 



1 



+ 



r(c)r(o + b - c) n _ ,„_„_,,_„_ 
r(a)r(6) " " • 



F\ [ c — a, 1 — a, c — a — 6 + 1; 1 



(|ar<7 z| < 7r, |ar<? (1 — z)\ < ir), 
is used. After this transformation, the hypergeometric function becomes 
2 Fi (2 - e IR , l + n x ,2 + n x , -£ w ) = {n x + l)\T(e IR - 1) 



T(n x + e /R ) 



, (1 + LY 1 *- 1 V 



+ 



Z=0 



£ 1 + f 



£ W J T{l + e IR )(n x -l)\ 



(3.3) 

Though the hypergeometric function in l.h.s. of Eq. (|3.3|) is finite when n x > 1 and ej R — > 0, 
the gamma function in r.h.s. has a 1/ejr pole. We have confirmed that the terms in brackets 
on the r.h.s. of Eq. (|3.3|) start 0(eir) when n x > 1. Thus, there is no 1/eir pole, as expected. 
Then, the remaining ^-integral in J4 becomes 

J4(s,t;0,0,0,0;n x ,n y ,n z ) = 

-1 / -t \ £IR 
^ t y y 4n ^2 J B ( n * + £ IR> n y + n z + e IR )n x \T(e IR )T(l - e IR ) 

1 / ~ \ —l—rix 



x / dw w n y- 2+£lR (i-w) r 



T(n x + ei R ) 



+ 



^ ( 1 



i=0 



The u>-integral can be solved in a term- by-term way for each power of I. The final form of J4 
is obtained to be 

J 4 (s,i;0,0,0,0; ——^—-B(n x + e ir , n y + n z + e IR )n x T{ei R )Y{l - eir) 



-t 



4TTflj 



(4tt) 2 s t 

: ' !f! -ty» />*;! • //-.//.,. • n y • sir) 
s 



T(n x + e ir ) 



x 2 F\ 1 + n x , n a + n y + e /jR , 1 + n x + n y + n x + e/ fl , 



4TTflj 



e IR n x 



i=0 



•1)' 



t J T{l + E IR ){n x -l)\ 



B(l + n y ,l + n z +£jr) 



x 2 F\ [1 + 1,1 + n z + e ir , l + l + n v + n z + e IR , 



(3.4) 



where t = t — iQ. 

When the numerator of the integrand is unity, = 0, the result of Ea. (|3,4j) is 

reduced to 

1 B(e ir ,£ IR )T(1 - £ IR ) 



J 4 {s,t; 0,0, 0,0; 0,0,0) 



(4tt) 2 s t 



£IR 



4tt^ 



sir 



2F1 [ l,em,l + em,~= ) + 



4^| 



sir 



U 

2F1 [ 1,£IR, 1 + EIR, — g 



(3.5) 
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We have checked that this result agrees with the precedence result obtained by Duplanzic and 
NizicjBJ in both physical and unphysical regions of kinematical variables, s and t. 

When n x = 0, the result of Ea. (|3.4|l has a shorter expression without a fake pole, as 



1 

(4vr~)2s t 



J<i(s,t; 0,0, 0,0; 0,n y ,n z ) = y—^- ^ B(e IR , n y + n z + e IR )T(l - sir) 



j B(l + n z ,n y + £i R ) 2F1 \l,n y + e IR ,l + n y + n z + e IR , r 



+ ^ 47r/ ^2 ) B(l + ny, n z + e ir ) 2 Fi ^1, n z + e/H, 1 + % + n z + e m , -| ) .(3.6) 

The infrared behavior of the loop integral can be obtained by expanding this formula with 
respect to £i R , as shown in Appendix B. 

In some cases with n x 7^ 0, we can avoid the fake pole by using the symmetry of the 
integrand. The basic integrand, Eq, (|3.2|) . is symmetric under the exchange x(n x ) <-* z(n z ). 
Then, the result with n x 7^ and n z = can be easily obtained as 

J4(s,t;0,0,0,0;n x ,n y ,0) = { J 4 (s, t, 0, 0, 0, 0; 0, n y , n z ), n z -» n x } 

for any values of n y . 

In the case of n x 7^ 0, n z 7^ 0, the integration has no IR-divergence with any value of n y > 0. 
However we cannot simply set si R to be zero, because there is the fake pole. After expanding 
Ea. (|3.4|l with respect to ei R , we confirmed that the ej R pole was canceled out. An explicit 
form of the Tyler expansion of the hyper geometric functions can be found in Appendix A. 

3.3 1 off-shell and 3 on-shell external legs 

When one of four external particles is off-shell, we can set it to be p\ without any loss of 
generality, p\ 0,p 2 =v\ = v\ = 0. Then, the integral of Eq. (|3.1|) can be 

Ji(s,t;pl,0,0,0;n x ,n y ,n z ) = 

T(2-e IR ) f 1 , f 1 -* , [ l - x -y x n *y n yz n * 
ax dy az- 



(4tt) 2 (4TTfJ, R Y IR J Jo Jo (-xzs - y(l - x - y - z)t - p(xy - iO) 2 £IR ' 

(3.7) 

After applying the same transformation as for the four-on-shell case, the integral becomes 



r(2 - em) 
(4tt)2 (4tt^) £ 



J±(s,t;pl,0,0,0;n x ,n y ,n z ) = — j-^B(n x + e IR , n y + n z + e IR ) 



-I ri w n y(l-w) n *v n * 
x dv dw- 



J 



(— s v(l — w) — t (1 — v)w — p\ vw — iO) 



The f-integral can be done as 



Z" 1 v n 
I v = dv 







t 2 {l + n 
where 



[(— s + s w + t w — p\ w) v — t w — iO] 2 £lR 
- 2-F1 ( 2 - e JR , 1 + n x , 2 + n x , , 



-i) " ir 



~ us 

t t W 
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Here, we use s + t + u = p\ . This result is the same as that in the four-on-shell case, except 
that u = — s — t + p\ instead of u = —s — t. After making the same transformation as that for 
the four-on-shell case, the remaining u>-integral can be expressed as 

J 4 (s,t;Pi, 0,0,0; 

) = ia \-? 7 B (. n x + £iR, n y + n z + e IR )n x \T(e IR )T(l - e IR ) 



x 



4vru| 



£IR 



(4vr) 2 s t 
T (i) 



s J T(n x + ei R ) 



4iru R 



v 



E 



(-1)' 



,(2) 



T(l + e IR ){n x -l)\ 



(3.8) 



where 



j(i) 

T(2) 



■/ 



w^+ny-^l* (1 - u,) n * ( 1 + 



s Jo 



— dw w riy (1 - w) Uz (1 + -W 



u 



u 

—w 

-l-l 



-l—n x 



t + u 
1 H z — w 



The first integration, jW, can be done as 

= 5(1 + n z ,n x + n y + e IR ) 2 Fi[l + n x ,n x + n y + e IR , l + n x + n y + n z + e IR , ~ 



u 
(3.9) 



For the second integration, 1^, we used a following binomial expansion: 



ki=0 



w 



Uy + kl 



t + u 

S 

Pi ~ s 



1 H w ) — 1 

s 



n v +ki 



E n y+kl C k2 (l + 

fc 2 =0 



t + u 



fe 2 



where m C n is combinatorial defined as 



m 



nl{m — ny. 

(2) 

Then, the second integration, X\ , can be done as 



(3.10) 



ki=0 
r-1 



pj-s 



E riy+k^ (-l) n v +k * ( — 
k 2 =0 V 



/ 

.70 



d-u; ( 1 + — iu 



-a+i) 



i + 



t + u 



-w 



k 2 +l-l+e IR 



n z riy+ki 

E E "s^l «y+felCfe2 (-l) fcl+fc2 
fei=0 fc 2 =0 



riy+ki 



1 



Z + fc 2 + £/R 



1 + T 



2 Fi M + /, I + k 2 + sir, 1 + 1 + k 2 + e IR , --= 

~ 2\ l+ k 2+em 
Pi 



Fi[ 1 + 1,1 + k 2 + e IR , l + l + k 2 + e IR , - 



ufi 2 
is 



(3.11) 
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where pi 2 = pf+iO. Here, the integral J^(s, t;p\ , 0, 0, 0; n x , n y , n z ) can be successfully expressed 
by a finite number of hyper geometric functions. The infrared structure of the tensor-integral 
can be obtained by expanding Eq. (|3.8|) with respect to Ei R . The results of expansions under 
the MS scheme are given in Appendix B. 

When the numerator of the integrand is unity, 
Eq. (j3.11|) is reduced to 



0, the second integration 



(In' [ 1 + — w 



2 Fi l,e IR ,l+e IR , 



1^ = z* 

s 
1 

SIR 



Then, the Eq. p.Hf) can be written as 



1 / i+u ^ £IR ~ 1 
1 + — —w 



pi 2 \ £ir p f. 1a . m 2 

2 F l l,e IR ,l + e IR , 



s J 



ts 



(3.12) 



Ms,t;pl 0,0,0; 0,0,0) 



1 B(e IR ,e IR )T(l - e IR ) 



(4vr) 2 s t 



Sl R 



£IR 



2 Fi [ l,e IR ,l+e IR ,-j ) + 



£IR 



u 

2F1 ( 1, ei R , 1 + ei R , — -z 



(-Pi 2 Y m 

J 4 (s,t; 0,0, 0,0; 0,0,0) 

1 B(si R ,e IR )T(l - e IR ) 
(4tt) 2 s t em 



2F1 I 1, ei R , 1 + ei R , 



upi 2 
is 



-ft 2 N e « 



2-^1 1, £TR, 1 + 



(3.13) 



We have again checked that this result agrees with the precedence result obtained by Duplanzic 
and NizicjS] in both the physical and unphysical regions of the kinematical variables. 

When n x = 0, the result of Eq. (|3.8|) does not have a fake pole. The result is obtained to be 



J 4 (s,t;p?,0,0,0;0,n y ,n 2 



T (l) 



+ 



(4vr) 2 s t 

\ &IR 



B(em, n y + n z + e IR )Y(l - e IR ) 



(3.14) 



where 
j(i) = 

t(2) _ 



B(l + n z ,n y + ei R ) 2 F x ( 1, % + £tr, l + n y + n z + e IR , 



n z riy+ki 
fe 1= fc 2 =0 



fci+fc 2 



i/ 



1 



&2 + S/J? 



Fi ( 1, & 2 + sir, l + k 2 + s IR , --=]-( — 



(3.15) 



2 Fi [l,k 2 + em, 1 + k 2 + sir, 



upi 



2\ 1 



ts 



(3.16) 



In the case of n x 7^ 0, n z ^ 0, the integration has no IR-divergence when using any value 
of n y > 0. However, we cannot simply set ei R to be zero again, because there is the fake pole. 
After expanding Eq. (|3.8j) with respect to ej R , we confirmed numerically that the sir pole was 
canceled out. The results of expansions under the MS scheme are shown in Appendix B. 
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4 Numerical check of the results 



The results of the vertex tensor-integral are rather trivial. However, those of the box tensor- 
integral are very complicated and highly non-trivial. We need some cross-checking of our results 
compared with other independent calculations. For the scalar integral of the one or two off-shell 
box integral, we can check our results numerically with those of the precedence calculation done 
by Duplanzic and NizicjH]. In both in the physical and unphysical regions of the kinematical 
variables, s and t, it was confirmed that the results given in this report agree completely with 
those in ref.[H]. 

The basic ingredients of a numerical calculation of the general case of the tensor integral 
are given in Appendix A. Those formulas are Eqs. (|A.3|) ~ (jA.5| ) to evaluate the hypergeometric 
function and Eas. (jA.10jW(jA.25j) to evaluate the generalized hypergeometric functions. 

At first, the results of a numerical evaluation of the hypergeometric function of the type 
Eq. ljA.4jl are compared with those obtained from Mathematical^ at several values of I, m, n 
and z at random. We confirmed that both results agree very well with each other to more 
than ten digits. For the generalized hypergeometric function, our recursion relation formulas, 
Eas. (|A. 10|) ~ (|A.25j) . were checked by comparing those of numerical integration of Ea. (jA.9j) 
using the numerical contour-integral (NCI) method jTlj developed by the author. The function 
F^(z), where n = 1,2 and 1 = 2, m = 3, given in Ea. (jA.9j) . was numerically evaluated at 
several values of z, as shown in Tables 1 and 2, which was compared with the NCI results. 
Both results gave very good agreement to about ten digits, as shown in tables. The imaginary 
part of the result must be zero, except z = 2000 + Oi case in the table. It was also numerically 
confirmed very precisely. 



z 


real/imag. 


analytic 


NCI 


2000 + 10~ i5 i 


real 
imag. 


3.167042847 x 10~ 7 
-1.308996938 x 10" 7 


3.167042848 x 10~ 7 
-1.308996939 x 10" 7 


-2000 + 10~ 15 i 


real 
imag. 


3.167042847 x 10~ 7 
C(10" 24 ) 


3.167042848 x 10~ 7 

e>(io- 18 ) 


0.2 + 10~ lb i 


real 
imag. 


3.723185361 x 10" 1 
e>(HT 16 ) 


3.723185362 x 10" 1 

e>(nr 12 ) 


-0.2 + 10 -15 i 


real 
imag. 


1.809694496 x 10" 1 
O(10" 16 ) 


1.809694496 x 10" 1 
O(10" 12 ) 



Table 1 Numerical comparison of ^3(2) between analytic formulas given in Eas. ljA.13^ ~ ljA.17fl and the nu- 
merical contour integral of Eq. lA.101 . 



z 


real/imag. 


analytic 


NCI 


2000 + 10~ 15 i 


real 
imag. 


1.024874615 x 10~ 6 
-9.949558053 x 10" 7 


1.024874616 x 10~ 6 
-9.949558053 x 10" 7 


-2000 + 10~ 15 i 


real 
imag. 


1.230491374 x 10 _ti 
C(10" 24 ) 


1.230491374 x 10" e 

e>(io- 17 ) 


0.2 + 10~ lb i 


real 
imag. 


1.624874690 x 10" 1 
e>(HT 15 ) 


1.624874690 x 10" 1 

o(io- 13 ) 


-0.2 + 10~ lb i 


real 
imag. 


1.005561296 x 10" 1 
C(10" 15 ) 


1.005561296 x 10" 1 
O(10" 12 ) 



Table 2 Numerical comparison of F^{z) between analytic formulas given in Eas. 1A. 211 ^ 1 A. 2511 and the nu- 
merical contour integral of Ea. ljA.210 . 
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n x fly n z 


real/imag. 


analytic 


NCI 


1 9 3 


V ' 1 1 

1 cell 

imag. 


—9 1 59Q8 x 1 0~ 9 
-2.78647 x 10" 9 


—9 1 5997 x 1 0~ 9 
-2.78650 x 10" 9 


2 2 


real 
imag. 


q 74570 x 10~ 9 
-3.22229 x 10~ 8 


q 74579 x 10 -9 
-3.22230 x 10~ 8 



Table 3 Numerical comparison of J,i(s, t; 0, 0, 0, 0; n x ,n y ,n z ) between analytic formulas given in Ea. ljB.55^ and 
the numerical contour integral. Here, we set kinematical variables in the physical region at s = 123, t = —200 
and /j,r = 1. 



t^x t^z 


real/imag. 


analytic 


NCI 


12 3 


real 
imag. 


-7.88683 x 10~ 1U 
-1.95176 x 10~ 9 


-7.88689 x 10~ 1U 
-1.95176 x 10~ 9 


2 2 


real 
imag. 


1.48133 x 10" s 
-2.04318 x 10~ 8 


1.48133 x 10" 8 
-2.04318 x 10~ 8 



Table 4 Numerical comparison of J4(s,t;pl, 0, 0, 0; n x , n v , n z ) between analytic formulas given in Ea. ljB.82fl 
and the numerical contour integral. Here, we set the kinematical variables in the physical region s = 123, 
t = -200 pi = 80 and p, R = 1. 

The numerical results of the box tensor-integral with zero and one off-shell external legs 
obtained using Eqs. (|B.55|) and (|B,82|) are also compared with those using the NCI method, 
integrating Eas. ()3.1|) and (|3,7|) directly. Both results show very good agreement within about 
five digits, as shown in Tables 3 and 4. 

5 Conclusions 

The general formulas of the 3- and 4-point tensor-integral were obtained directly without any 
reduction method to the scalar integrals. The IR behavior of the tensor integrals was clearly 
shown by expanding the results with respect to the extra space-time dimension due to the 
dimensional regularization. All results were expressed by very short formulas in a suitable 
manner for a numerical calculation. The results of the scalar-integral were compared with the 
precedence results, and showed complete agreement in both physical and unphysical regions of 
the kinematical variables. For the IR finite case, the analytic results were compared with the 
numerical contour integration, and gave a consistent result within the numerical integration 
error. 
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Appendix 



A Numerical calculation of the hypergeometric function 

In this Appendix, the basic properties of the hypergeometric function^] and their numerical 
evaluation are summarized. The 1/eir expansion of the hypergeometric functions appearing 
in the tensor-integrals and their numerical evaluation are also shown. 



A.l Hypergeometric function 

The Gauss-series representation of the hypergeometric function is 

2 F 1 (a,b,c;z) = 2 Fi(b,a,c;z) 



fc=0 v ; 



where (•)& is Pochhammer's symbol defined as 

Fja + k) 
{a)k = ~T(a)~- 

The Euler integral representation is 

2 F 1 (a,b,c;z) = ^ ( . c) C T b ~\l - r) M (l - ZT)~ a dr, (A.2) 



r(6)r(c-6) J 

(Kc >TZb> 0). 

When a is a negative integer, such as a = —m, the Gauss series is truncated at k = m, and 
becomes a polynomial, 



iFi(-m,b,c;z) = ^ 



(-m) k (b) k z k 
(c) k k\ 

k=0 y ' 



= jt^M-zf, (A.3) 
fc=o {c)k 

where m Cfc is combinatorial defined in Eq. (|3,l()|) , For numerical evaluations of tensor- integrals, 
the following type of the hypergeometric function might be numerically calculated as 

2 Fi(l,m + l,n + m + 2; z) 

Wi)i V(l - r) " (1 - zT) ~' iT 

y y ^i+k, a m +kl c k i r {1 _ ZT) -i+k 2d (A 4) 

^ > B(m + l,n + i) z m+k ^ Jo 

fc 1= k 2 =0 J[) 

where I, m, n are positive integers. Here, integration can be performed as 

i r_Miz£) k2 -i + i = o 



B(m + 

n m+ki 



{l-zT)- l+ ^dT= I . 2 fl , fc2 _ i+1 1 (A.5) 
This formulas can be used for a numerical evaluation of hypergeometric functions of this type. 
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A. 2 Generalized hypergeometric function 

For a Laurant expansion of the hypergeometric function with respect to sir, the following 
generalized hypergeometric function is necessary: 



3F 2 ({a 1 , a 2 ,a 3 }, {h, b 2 }; z 

n * 

(bi)n(b 2 ) n n\ 



\fl\ )n 

2-. (h,\(h^ ^ A - b > 



n=0 

r(6i)r(62) 



r(oi)r(6i - oi)r(o2)r(62 - 02) 

x jf (It J dv v ai - b2 (l-v) bl - ai - 1 T a2 ~ 1 (v -T) b2 ^ a2 - 1 (l- zt)' 113 , (A.7) 



and 



4,F3{{a 1 ,a 2 ,a3,a i },{biMMY,z) = > f . (A.8) 

^ (Oljn(02jn(,03jn 

Those generalized hypergeometric functions appear in a following integral: 

P hU z ) = ^ff Q dT r jl -\l-ZT)-& + V In™ r, (A.9) 

where ji is a positive integer and j'2 is an integer (it can be negative). 
When n = 1, the integral becomes 



^■ 2 (*) = " £ dr r^-\l - zr)-^ +1 ) lnr, 



'0 

3-^2 ji, h + + 1, ji + 1};^) 



(A.IO) 
(A.ll) 



5 (ji + mn + *: + i)so"2 + 1, * + 1) ■ (A,12) 

When j'2 is a negative integer, such as j 2 = —j < — 1, this function can be express by a 
polynomial, 

When Ji = 1 and j'2 = 0, 1, the integral can be performed easily as 

*8!m = (A.i4) 

F<V W = -!^. (A.15) 
When ji = 1 and j 2 > 1, we can use the following recursion relation: 

^ _ h m , (i-^)- i2 -i 



For the general case, the function F^- 2 (z) can be obtained using F^ 1 ) , 



= ^D-l^-iC*^). (A.17) 



fc=0 
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When n = 2, the integral becomes 

X - ^ dr T jl -\l - zr)-( J2+1 ) ln 2 r 
4F 3 ({j 1 J 1 ,j 1 J 2 + 1}, {ji + 1, h + 1, ji + i}; z) 



F {2) . (z) 



E 

fc=0 



(ii + fe) 3 (j 2 + fe + l)B(j 2 + l,A;+l) 



When j2 = — j < — 1, it is also represented by a polynomial, 



= e 



fc=0 



i-iC k (~z) k 



When J\ = 1 and j'2 = 0, 1, the integral can be performed easily as follows: 

f.(2)/^ _ Li s(^) 



z 

Li 2 (z) 



When ji = 1 and j 2 > 1, we can use the following recursion relation: 



02 + i)^ +1 W-M^W-f^W 



£■(2) 



0. 



',(2) 



For the general case, the function j- 2 (z) can be obtained using F{/(z), 



ji-i 



(A.18) 
(A.19) 

(A.20) 

(A.21) 

(A.22) 
(A.23) 

(A.24) 
(A.25) 



k=0 



A. 3 expansion of a hypergeometric function 

In the general form of the tensor integrals, following type of integral appears: 

1l,m,n = f r l+n - 1+ ^(l-Tr(l-ZT)-( l+1 Ur, 

Jo 

= B(l + mJ + n + s m ) 2 Fi(l + l,l + n + e IR , 1 + I + n + m + e IR , z) . (A.26) 

In order to show the IR structure of the function, a series expansion of the function T[ m n with 
respect to sir might be considered, 



(A.27) 



where l,m,n are non-negative integers. When / = m = n = 0, the Laurant expansion of the 
function can be 



T r (1 - 



1 (1-zt)- 1 



-sir 



dr, 



B(l, Sir) 2F1 (1, Sir, 1 + Sir, z) 



f j ^ x-i {l-zry 1 -I lnr /. \ 
= / ^ ) l + - + eiR Ul-zr -1 



_±_ _ i n (i _ z ) _ Li(z) £/i? + O ( £ y . 
sir 



(A.28) 
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Then the first three terms of Eq. (jA.27j) are 



•^0,0,0 ( z ) ~ 1j 
4%o(z) = -ln(l-^) 
4o,oW = -U(z). 



When l = ra = 0,n^0, the integral becomes 

X 0A n = /V-^a-zr)- 1 ^. 
Jo 

Then, the first three terms of the Eq, (|A.27|) are 



T (o) / \ _ 2Fi(l,n,l + n;z) 

J 'Q,0,n\ z ) — 

•^0,0,n( Z ) = — ^n,d( Z )' 



When I = 0, m ^ 0, n = 0, the integral becomes 



X0,m,0 = / T- 1+£l *(l-T) m (l-ZT)- l dT 
JO 

m .x 
= J2 m C k (-l) k T k - 1+S ™{l-ZT)- l dT 

7- „ JO 



k=0 
in 



— E rnCk( — l) k 1o,0,k- 
k=0 

Then, the first three terms of the Ea. (|A.27|) are 

m 

F 3,m,o( Z ) = E "iC fc (-l) fc JF 1(j 



k=0 

m 

fe=0 

When l = 0,m^0,n^0, the integral becomes 



*gU*) = E mft(-l)*4(z). 



'-0,m,n 



I T »-l+=«(l - T ) m (l - zr)- 1 ^ 
V m C fc (-l) fc / r^-^Cl-^r)- 1 ^ 

m,Ck{ — l) fc 2o,0,n+fc- 



fc=0 

Then, the first three terms of the Eq. (jA.27|) are 

^§!»,n(*) = B(l + m,n) 2^1(1,71, 1 + n + m;^) 

fc=0 
m 

fc=0 
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When I ^ 0, m = 0, n = 0, the integral becomes 

3,0,0 = C T l - 1+eiR {l-ZT)-^dT. 

Jo 

Then, the first three terms of the Ea. (|A.27|) are 



( Z ) - t 1 '^ 1 
v {\) i , 2 Fi{l,l,l + l;z) 



When Z ^ 0, m = 0, n ^ 0, the integral becomes 

l,n n = [ T l+n - 1+£ 'K(l-ZT)-V + VdT. 







Then, the first three terms of the Ea. (|A.27|) are 



T (o) , s 2-Pi(l + /, I + n, 1 + I + n; z) 



l + n 

= 

When / ^ 0, m ^ 0, n = 0, the integral becomes 

2i,m,o = / r'- 1+£ ^(l-rr(l-^rr^ +1 )dr. 

JO 

Then, the first three terms of the Ea. (|A.27|) are 

^S,o(*) = + H 2^1(1 + + / + m; 2 ), 

fc=0 
m 

*#U*) = E m^(-i)*^w. 

fc=0 

When l/0,m^0,n^0, the integral becomes 



k,m,n = j 1 T l+n - 1+£lR (l-T) m (l-ZT)-^dT. 







Then, the first three terms of the Eq. ([A.27|) are 

tfmn = B(l + n,l + m) 2 F 1 (l + l,l + n,l + l + m + n;z), 



" r l,m,n ~ 2—i mUfe v L ' r l+n+k,l\^h 
k=0 

m 

• r l,m,n ~ 2—/ 1 ) r l+n+k,l\^ ] 

k=0 



(A.48) 
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B 1/eiR expansion of tensor-integrations with the MS scheme 



We give 1 /sir expansions of tensor-integrals under the MS scheme in this Appendix. The MS 
scheme is realized with the following replacement of the renormalization energy-scale: 



An 



where 7^ is Euler's constant. 

B.l Vertex with a 1 off-shell external line 

The general result of the vertex tensor-integral is shown in Eq. (|2.2|) . The IR structure of this 
function can be obtained by a Laurant expansion of the beta function for the IR singular case. 



B.l.l 



71,, 







When n x = n v = 0, the tensor integral becomes 



J 3 (0,0,p|;0,0) 



r(-ejfl) 
(An) 2 



-P3 



2 \ £ IR 



Anp 



After the ejr expansion with the MS scheme, 



J 3 (0,0,p§;0,0) 



1 



(An) 2 p{ 



C 



(-2) 
31 



+ 



'IR 



(-1) 
31 

£IR 



1 B(eir,Eir) 

-pi 2 



+ C^+0(e IR ) 



where 



r (-2) 


= 1, 


°31 


= ln(- 


r (o) 

u 31 


n 2 
12 



A 



+ - In 2 



-vl 

/4 



(B.l) 

(B.2) 

(B.3) 
(B.4) 

(B.5) 



B.l. 2 n x + n y ^0 

The general result of the vertex tensor-integral, Ea. H2.2|) . is symmetric under n x 
changes. When one of n x or n y is non-zero, the tensor integral becomes 

J 3 (0,0,p 2 ;n,0) = J 3 (0,0,p 2 ;0,n) 
£irT(— Eir) 



and n y ex- 



(An) 



-p 3 2 Y IR 1 B(n + E IR ,E IR ) 



Anp\ 



-pi n + 2e ir 



(B.6) 



After an eir expansion with the MS scheme, 

J 3 (0,0,^71,0) = J 3 (0,0,p§;0,n) 



1 



(An) 2 p 2 



C 



(-1) 

32 



£IR 



+ cj$ + 0(e IR ) 



where 



C 



(-1) 

32 

(0) 
32 



1 

n 



2 1 , - 
-- 2 + - ( In 



A 



(B.7) 

(B.8) 
(B.9) 



where H m is the harmonic number, defined as 

-u v-^m 1 

rI -m — 2_f7=l 7 > 
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B.1.3 n x ^O,n y ^O 

In this case, because there is no IR-divergence in Ea. (|'2.2l) . em can be set to zero, 
^3(0,0,^;^,^) 



1 (n x - l)\(n y - 1)! 



(B.10) 



(4vr) 2 p 2 (n x + n y )\ 

B.2 Vertex with two off-shell external lines 
B.2.1 n x = riy = 

The function Go(z), defined by Eq. (j2,4j) . can be expanded with respect to em around zero as 

g (z) = G^{z)e IR + Q ( i\z)em + 0{em), (B.ll) 

where 

ln(l - z) 



ln 2 (l - z) 
2z 



(B.12) 
(B.13) 



Then, the vertex integral Js(0,P2,pi',0,0) can be express after an em expansion with the MS 
scheme as 



J 3 (0,p|,^;0,0) 



l 



(4vr) 2 pi 



C 33 , AO) 



£IR 



+ C^+0(e IR ) 



where 



(-1) 



»(-2)-(l) (P3-P2 



33 



AO) 
u 33 



u 31 ^0 



L 31 ^0 \ ~2 



~ 2 / ' 
P3 

.2 ^2 



(-i)rr(i) (P3-P2\ ,A-2)rP) Pzzll 



-t- C-3! f 



,2 „2 



P3 2 



(B.14) 

(B.15) 
(B.16) 



B.2. 2 ^ = 0,71^/0 

Because of vertex integral Js(0,P2,p^) 0, n) is IR finite, we can set em to be zero as 

</3(o,P2>p!;°> n ) 

where 



1 AO) 

2 „2 U 34 i 



(4vr) 2 pi 



a 



(0) 

34 



r (-l) c (l) / gj ~ Pi 
u 32 ^0 ~ 2 

Pz 



(B.17) 



(B.18) 



B.2. 3 n x ^0,n y = 

For non-zero values of n, Q n {z) can be expanded with respect to e/^ around zero as 

g n (z) = g(?)( z )+g(i)( z )e IR + G(e 2 IR ), (B.19) 



where 



G^U 



n 



n + l 
1 

n + l 



2^(1,1,2 + 7^), 

2 Fi(l,l,2 + n;z) -n 



z-iy in 2 (l-z) 
2z 



k=l 



ra _ fc (l-z) fc (l-fcln(l-z))-l 
/c 2 z 



(B.20) 



(B.21) 
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Then, the vertex integral Js(0,P2,p^; n, 0) can be express after an sir expansion with the MS 
scheme as 



^3(0,p 2 ,p 3 ;n,0) 



1 



(4vr)2 pi 



C 



(-1) 

35 
£IR 



+ 4° 5 ) + 0(e IR ) 



(B.22) 



where 



C 



(-i) 

35 



(-i) c (o) ( Pi-Pi 

32 ^0 \ p ^2 



AO) 
u 35 



32 »0 



(-i)rj(i) r gj-Pg 



P3 



(o)-(o) / v\-v\ 



' 31 »0 



P3 2 



(B.23) 
(B.24) 



B.2.4 n^O.n^O 

In this case, there is no IR-divergence and ejr can be set to zero, 

1 (n x - l)l(n y - 1)! 



J3(0,pl,pl;n x ,n y ) 



(n x + n y )\ 



2 2 

x 2^(1,1,2 + ^;^-^ 

P3 



n, 



n x + 1 



(B.25) 



B.3 Box integral with all on-shell external legs 

The general result of the box tensor- integral with all on-shell external legs is shown in Eq. (|3,4j) . 
The IR structure of this function can be obtained by a Laurant expansion of the beta function 
for the IR singular case. 

B.3.1 n x = riy = n z = 

When the numerator of the integrand is unity, n x = n y = n z = 0, the box tensor-integral 
result, Ea. (|3.4[) . is reduced to 



J 4 (s,t, 0,0, 0,0; 0,0,0) 



1 B(ei R ,e IR )T(l - e IR ) 



(4vr) 2 s t 



sir 



47T/U 



s \ £IR ( u\ ( -t x £IR 



Fi l,e IR ,l + ei R , 



ii 



(B.26) 



After using a Laurent expansion of the hyper geometric function of Ea. (|A.28|) the loop integral 



with the MS scheme is 



Ji(s,t, 0,0, 0,0; 0,0,0) 



C 



C 
C 



(-2) 
41 

(-1) 



41 

(0) 
41 



1 



(4vr) 2 s t 
4, 
2 

J2 



r {-2) r (-l) 

41 +^+C^+0(e IR ) 



s 



IR 



£IR 



(B.27) 

(B.28) 
(B.29) 

-T- 2L H-ih 2L H-Tj+n^) +ln2 (^ 

-21n ( -j) In (l + I) -2 In In ( 1 + % ). (B.30) 



In I -£■) +ln / 



7T 
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B.3.2 n x = 0,n y = 0,n z ^O 

When n x = n y = 0, n z / 0, the result is reduced to 

B(eir, n z + e IR )T(l - e IR ) 



J 4 (s,t; 0,0,0,0;0,0,n z ) = 



(4tt) 2 s t 



sir 



B(l, n z + ei R ) 2 F\ 1, n z + e/ R , 1 + n z + e /jR , - 



+ 



(B.31) 



Laurent expansions of beta- and gamma-function with the MS scheme are obtained as follows: 



B{ni + e IR , n 2 + e IR )T(l - e IR ) 



2 \ - IB A-l) / 2\ 

^ l?, +4n 1 (? 2 )+^,«(9 2 )^ 



where 



= In 



4V) 



/4 



and 



•^11,712(5 ) 

•^11,112 ) 



3W£_ 1 + 7r 2 + 3 In 



0. 



In 



2H n -9^ {1) ( 



7) 



B(m,n 2 

B(ni,n 2 ) ( H ni -i + W n2 - 2H ni+ „ 2 _i + In 



>4 



(B.32) 

(B.33) 
(B.34) 

(B.35) 

(B.36) 
(B.37) 

(B.38) 



where n, n\, n 2 are positive integers. Here, ip^(z) is the first derivative of the digamma 
function, given by 

2 n ~ 1 1 
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Finally, the tensor integral can be obtained to be 
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B.3.3 n x = 0,n y ^ 0,n z = 

When n x = and both and n 2 are non zero, J4 becomes 
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One can see that Ji(s, t; 0, 0, 0, 0; 0, n y , n z ) is symmetric under a simultaneous exchange of 
t <-> s, n y «-» ra 2 . Then, for the case of n x = n 2 = 0, 7^ 0, the required integration can be 
obtained as 

J 4 (s,i; 0,0,0; 0,^,0) = { J 4 (s, t; 0, 0, 0; 0, 0, n z ), n z -> n y , s «-► t}. (B.44) 



B.3.4 n x = 0,ra y / 0,n z 7^ 

Let's start from ea. (|B.43|) . There is no IR-singular pole in the square-bracket and the ^ pole 
in the Beta function in front of the square brackets. It is thus needed to expand the terms in 
the square-bracket up to O (em) . 
Then, finally we obtain 
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After an em expansion with the MS scheme, it is obtained as 
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B.3.5 7^ 0, n y = 0, n 2 = 



The basic integrand, Eq. (j3.2j) . is symmetric under the exchange x(n x ) <-» z(n z ). Then, the 
result can be easily obtained as 



J 4 M;0,0,0,0;n x .,0,0) = { J 4 (s, t; 0, 0, 0, 0; 0, 0, n 2 ), n z -► n x }. 
B.3.6 n x ^0,n y ^0,n z = 

When n 2 = 0, the result of the integration, Ea. (|3.4[) . can be 
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After an £/# expansion with the MS scheme, it is obtained as 
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B.3.7 n x / 0, n v =&ny,n z ^ 



In the case of n x ^ 0, n z ^ 0, the integration has no IR-divergence with any value of n y > 0. 
The result is given by 
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B.4 Box integral with 1 off-shell and 3 on-shell external legs 



The general result of the box tensor-integral with 1 off-shell and 3 on-shell external legs is 
shown in Eq. ()3.8|) with Eas. (|3,9|) and ()3.11|) . The IR structure of this function can be obtained 
by the Laurant expansion of the beta function for the IR singular case. 



B.4.1 n x = n y = n z = 

When the numerator of the integrand is unity, = 0, the second integral call^ is 
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where p\ = p\ + iO. Then, the final result can be written as 
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After an e/^j expansion with the MS 1 scheme, it is obtained as 

1 



J4(s,t,p{, 0,0,0; 0,0,0) 



(4tt) 2 s i 



r (-2) ^(-1) 

^ + ^+cJS ) +0(e«) 



sir 



where 



.(-2) 
"46 

.(-1) 
"46 



^(0) 
--46 



7 (-2) 
"41 

c (-i) _ 2 

L 41 Z 



In 



In 1 + 



Is 



c (o) , 

L 41 + 



7T- 



ttj) 2 



— +2Li 2 ( — ) _ln 



2 r "Pi 

/4 



+ 21n( 



-Pt 
V Mr 



iu!i + | 



(B.59) 

(B.60) 

(B.61) 
(B.62) 

(B.63) 



B.4.2 n x = 0,n y + n z > 1 

When n x = 0, and n y , n z have any value with n y + n z > 1, the final result can be written as 
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After an £/# expansion with the MS scheme, it is obtained as 
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B.4.3 n x ^ 0, n y = 0, n z = 

When n x 7^ 0, and n y = n z = 0, the final result can be written as 
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After an £j R expansion with the MS scheme, it is obtained as 
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B.4.4 n x ^ 0,n y + n z > 1 

At last we treat the most general case with n x ^ 0, and n y + n z > 1. The final result can be 
written as 
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In this case, the result might be IR finite. Then, after an ei R expansion with the MS scheme, 
it becomes 
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